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一 、 


要 ;地面 颤 振 试验 (ground flutter test,GFT) 系统 采用 少量 激 振 器 模拟 连续 分 布 气动 力 ,为 实现 


气动 插值 点 的 缩聚 ,推导 了 地 面 颤 振 试验 系 统 的 控制 方程 ,对 比 了 GFT 系统 和 原 颤 振 系 统 控制 方程 
在 广义 力 项 区 别 ,定义 了 激 / 测 振 点 位 置 优化 的 目标 函数 ;为 了 提高 插值 精度 ,充分 利用 结构 振动 已 
知 信息 ,提出 了 虚拟 激 / 测 振 点 技术 ;采用 分 群 粒 子 群 算法 兼顾 局 部 寻 优 和 全 局 寻 优 ,搭建 优化 流 


程 ， 


对 插值 点 进行 优化 配置 。 基 于 平板 机 机 开 展 了 GFT 系统 激 / 测 振 点 位 置 优化 计算 ,构建 了 地 面 


集中 气动 力 并 进行 了 颤 振 特 性 测试 试验 。 试 验 结果 表明 ,本 研究 提出 的 方法 精度 较 高 ,满足 地 面 颤 


振 试 验 的 需求 。 
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Sensor and shaker locations optimization 
of the ground flutter test system 
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;Abstract :The ground flutter test (GFT) system is the testing equipment which simulates the continuously 


distributed aerodynamic force with a small number of shakers. To realize the polycondensation of aerody- 


namic-force interpolating point, the GF'T’s control equation is derived and the generalized-aerodynamic- 


force difference between the GFT system and the original flutter system is showed. Based on the equiva- 


lence of generalized aerodynamic force, the objective function about sensor and shaker locations optimiza- 


tion is defined. The method proposed to improve the interpolation accuracy of the mode of vibration is 


named virtual sensor-shaker technology. Multi-Swarm Particle Swarm Optimization Algorithm ( MSPSO) is 


used to establish the optimization process,which has a good performance to moderate the contradiction of 


global searching ability and local searching ability. The GFT about a plate wing case was carried out, in 


which the sensor and shaker locations were optimized with the method proposed in this paper. The test re- 


sults verified that the accuracy of the optimization method is high enough for the GFT. 
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飞行 器 地 面 颤 振 试 验 ( ground flutter test, GFT) 
技术 是 一 项 利用 集中 力 模 拟 连续 分 布 气 动力 ,在 地 
面 进 行 颜 振 测 试 的 全 新 前 沿 技术 。 该 技术 在 地 面 振 
动 试验 基础 上 发 展 而 来 ,其 原理 是 通过 布置 的 传 感 
器 (加 速度 .速度 和 位 移 ) 测试 结构 的 振动 响应 , 通 


如 速 此 ， 和 这 税 信 品 


忆 改 项 树 


功率 虐 六 监 


过 提前 重 构 的 气动 力 模型 进行 气动 力 实时 计算 ,并 
利用 激 振 咒 将 等 效 气动 力 加 载 至 结构 上 ,从 而 模拟 
结构 的 颤 振 特 性 ,实现 颤 振 边 界 的 测试 和 验证 ,具体 
如 图 1 所 示 … 。 


实时 浊 定 党 气动 力 MIMO 


洲 助 思 


控制 营 


控制 出 扑 


| 


dsPACE 


习 1 地 面额 振 试 验 技术 原理 图 


Fig.1 Diagram of the ground flutter test system 


全 早 在 20 世纪 60 年 代 , Kearns'2 对 地 面 颜 振 模 
报 谤 行 了 初步 研究 。2001 年 ,俄罗斯 TsAGI 提出 了 
EMKM( 电动 机 械 模拟 ) 技术 ,计算 机 根据 实时 测 得 的 
结 梅 振动 信息 实时 计算 气动 力 ,并 通过 激 振 器 加 载 
到 结构 上 , 实现 了 水 平 尾 经 的 地 面额 振 模 拟 中 。 
20t 年 ,美国 ZONA 公司 Zeng 等 “提出 了 干 风 洞 
的 申 想 ,文中 对 激 / 测 振 点 位 置 进行 了 优化 配置 ,并 
引 愉 MIMO 控制 器 对 力 加 载 进行 控制 ,实现 了 对 和 矩 
形 银 板 机 翼 的 地 面 颜 振 模拟 。20 世纪 80 年 代 , 国 
内 的 潘 树 祥 等 "9 开展 了 用 地 面 共 振 方法 进行 热 颜 振 
试验 的 研究 ,当时 采用 的 是 模拟 计算 机 控制 ,对 超 音 
速 非 定常 气动 力 的 模拟 采用 活塞 理论 计算 ,并 直接 
简化 到 机 翼 的 6 个 气动 分 区 ,用 6 个 激 振 器 来 施加 
模拟 的 气动 力 。 北 京 航空 航天 大 学 许 云 涛 等 ”开展 
了 地 面额 振 试验 中 的 非 定常 气动 力 模拟 研究 ,文中 
采用 广义 力 相 等 的 方法 对 缩减 的 气动 力 加 载 位 置 进 
行 优化 ,并 进行 了 颜 振 时 域 仿真 计算 。 

地 面额 振 试 验 的 核心 关键 是 获得 气动 插值 点 的 
位 置 ,并 在 此 基础 上 建立 非 定常 气动 力 的 重 构 模 
型 中 。 气 动 插值 点 包括 测 振 点 和 激 振 点 ,两 者 分 别 
对 应 着 颤 振 分 析 对 象 的 连续 弹性 变形 离散 缩减 以 及 
连续 分 布 气动 力 的 离散 缩减 ,该 过 程 会 引起 弹性 变 
形 测量 和 气动 力 加 载 的 原始 误差 ,因此 为 了 在 少量 


的 激 振 点 和 测 振 点 数量 下 控制 缩聚 的 误差 ,需要 对 
激 ( 测 ) 振 点 数量 及 位 置 进行 优化 ,符合 误差 要 求 的 
激 ( 测 ) 振 点 是 地 面 闸 振 试验 的 后 续 实 施 前 提 '” 。 
当前 , 激 振 点 和 测 振 点 优化 目标 定义 方式 主要 
有 两 种 ;一 种 是 直接 以 颤 振 速度 和 颤 振 频 率 为 
优化 目标 ; 另 一 种 以 广义 气动 力 相 等 为 优化 目标 。 
第 一 种 方法 相 比 第 二 种 方法 ,目标 直接 明了 ,但 是 计 
算 量 较 大 , 且 由 于 缺少 对 系统 进行 内 在 的 机 理 分 析 ， 
系统 赢 振 分 支 可 能 发 生变 化 ,同时 与 地 面 阁 振 试验 
测试 额 振 边 界 的 目标 相悖 ; 以 广义 气 劲 力 为 优化 目 
标的 方法 ,优化 计算 量 小 ,需要 对 系统 的 颤 振 机 理 进 
行 分 析 , 选 取 振 型 为 优化 目标 ,并 依据 赢 振 参与 量 对 
振 型 进行 加 权 处 理 ” ,能 够 保证 闸 振 机 理 的 一 致 性 。 
本 研究 通过 推导 GFT 系统 的 闸 振 方程 ,前 明了 基于 
广义 力 相等 的 激 / 测 振 点 优化 技术 理论 原理 ,提出 了 
虚拟 激 / 测 振 点 技术 ,减少 了 激 / 测 振 点 个 数 , 提 高 了 
气动 力 缩减 的 精度 ,并 基于 机 翼 模 型 开展 了 气动 力 
模型 仿真 与 试验 测试 ,证 明了 本 研究 方法 的 有 效 性 。 


1 地 面 颤 振 系统 控制 方程 推导 


测试 对 象 气动 外 形 发 生变 形 时 ,CFT 系统 受到 
结构 惯性 力 、 结 构 弹 性 力 、 激 振 絮 实时 激 振 力作 用 ， 


在 有 限 元 节点 自由 度 集 x 下 的 时 域 颤 振 控制 方程 可 
以 写 为 2.9 


Mx +Kx -f=0 (1) 
式 中 :MK 分 别 为 质量 矩阵 .刚度 矩阵 ;xz 为 有 限 元 
节点 自由 度 集 ; 拓 为 与 上 等 效 的 x 自由 度 集 下 的 力 。 
由 虚 功 相等 原理 可 得 关系 式 记 = Gf., 代 入 式 (1) 
可 得 


Mx +Kx -Gf.=0 (2) 
式 中 :f 为 激 振 点 自由 度 集 e 下 的 激 振 力 ;G 为 x 
自由 度 集 到 e 自由 度 集 的 样 条 插值 矩阵 。 假 定 空气 
动力 网 格 中 心 点 的 自由 度 集 为 4a, 对 应 的 气动 力 癌 
量 为 ,缩减 得 到 激 振 力 ,代入 式 (2) 得 
MX +Kx -GG,.f,=0 (3) 
基于 小 扰动 假设 ,时 域 气动 力 f 的 表达 式 为 


£1) = Je-[a(7 -asdar (4) 
式 虽 :如 为 气动 力 传递 函数 ;9。 为 无 限 来 流动 压 ;、 
VO@b 为 参考 长 度 ( 取 参 考 弦 长 的 一 半 ) 和 无 限 来 流 
束 度 ;代入 式 (3) 并 作 拉 普 拉 斯 变换 得 
CMx (5) + Kr(s) -9。 Gr6, [al Ja) -0 
NN (5) 
四 引入 模 态 坐 标 4 和 结构 模 态 振 型 矩阵 B,、 罗 ,， 
帅 模 态 释 加 原理 有 x = 办 da = p,q, 代入 式 (5) 并 左 
乘 友 得 
‘M+R G0 [BY)]b. hs) =0 
- (9) 
式 钙 :M 和 天 分 别 为 广义 质量 矩阵 和 广义 刚度 所 
阵 。 引 入 关系 式 办 = Gg,,G,, = G.G,。, 式 (6) 继 
续 可 以 写 为 

(SM+R-g.0 G76. [BY ) sds) =0 

a 

上 式 即 为 假设 结构 变形 准确 测 得 ,并 只 考虑 由 
分 布 气动 力 缩减 到 有 限 个 激 振 力 情况 下 的 GFT 系 
统 颜 振 控制 方程 。 实 际 上 ,基于 有 限 个 测 点 测 得 的 
变形 也 是 存在 误差 的 , 测 点 位 置 构成 自由 度 集 变 
形 为 c, 则 由 测 点 测 得 的 带 误差 的 变形 为 w' = Gue， 
c = Gua, 代 入 式 (7) 可 得 考虑 结构 变形 插值 误差 的 
完整 CFT 系统 额 振 控制 方程 , 即 


(S21+R -9.6 G16. [AB(Y)]6.6.¢,)a(s) =0 
(8) 
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为 了 和 颤 振 系 统 的 颜 振 方程 进行 比较 ,给 出 闫 
振 系 统 的 颤 振 方 程 , 即 
(SM+k-gp [a(Y)]g.)a(s) =0 (9) 
可 见 两 者 的 差别 在 于 广义 气动 力 项 上 。 相 较 于 
颤 振 系统 ,DWT 系统 的 颤 振 方 程 的 空气 动力 转换 矩 
阵 五 (sLAV) 分别 左 乘 了 (GG,,)"、 右 乘 了 GG,。 
这 两 项 分 别 对 应 力 减 缩 和 振 型 扩充 引入 的 ,通过 激 
振 点 和 测 振 点 位 置 优化 ,让 两 者 广义 气动 力 差 别 减 
小 ,就 是 基于 广义 气动 力 相等 的 激 / 测 振 点 位 置 优化 
技术 。 


2 优化 目标 函数 定义 


由 式 (7) 可知, 在 不 考虑 变形 测试 误差 的 情况 
下 ,为 了 使 广义 气动 力 相等 , 则 有 
b=G.G..0, (10) 

上 式 的 物理 意义 是 :将 激 振 点 自由 度 表示 的 振 
型 插值 到 气动 网 格 节点 自由 度 时 ,得 到 的 振 型 和 原 
气动 网 格 节点 自由 度 下 的 振 型 一 致 , 即 激 振 点 构成 
的 自由 度 可 以 准确 地 表示 原 系统 的 振 型 。 

如 果 激 振 点 构成 的 自由 度 能 够 准确 表示 某 一 阶 
振 型 , 则 GFT 系统 模拟 的 该 阶 振 型 广义 力 是 准确 
的 ,各 阶 模 态 振 型 误差 对 广义 力 的 影响 是 相互 解 耦 
的 。 对 于 赢 振 系统 而 言 , 对 疗 振 贡献 度 较 大 的 主要 
为 低 阶 模 态 ,因此 只 要 插值 点 能 够 准确 表示 低 阶 模 
态 的 阵型 , 则 得 到 的 颜 振 边界 就 与 原始 颜 振 边界 接 
近 。 为 了 对 激 振 点 分 布 进行 优化 ,首先 需要 对 原 系 
统 进行 传统 的 频 域内 颤 振 分 析 , 确 定 系统 的 第 一 癌 
振 分 支 的 主 耦合 振 型 ,选取 主 耦合 振 型 以 及 其 它 贡 
献 量 较 大 的 振 型 作为 优化 目标 振 型 ,形成 振 型 矩阵 
4 ,计算 如 下 目标 函数 。 

_ 168. -81 _ 16- 
"gl 区 沁 | 


(11) 


EY 4 -bo 


fo ,min 
/这 六 同 
目标 振 型 加 权 , 由 于 选取 的 目标 振 型 中 , 主 耦合 
振 型 的 插值 精度 是 需要 优先 保证 的 ,因此 需要 对 所 
选取 的 目标 振 型 进行 加 权 处 理 ,重要 的 振 型 分 配 大 
的 权重 系数 。 加 权 后 的 目标 函数 为 


(12) 
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到 的 气动 网 格 自由 度 振 型 。 为 了 使 GFT 系统 的 广 


EE tb) 


535 
测 振 点 的 目的 是 准确 测 出 结构 的 变形 , 由 式 
(8) 可 知 ,为 了 使 广义 力 相 等 , 则 有 
p= CC (14) 
测 振 点 总 的 优化 原则 同样 是 使 测 振 点 的 位 置 能 
准确 地 表示 系统 振 型 ,因此 和 激 振 点 的 优化 目标 函 
数 形式 一 样 ,这 里 不 再 蒙 述 。 


op ( 13) 


3 ”虚拟 激 / 测 振 点 技术 


对 于 具有 已 知 边界 条 件 的 模型 ,为 了 提高 插值 
精 苦 ,采取 通过 增加 虚拟 激 ( 测 ) 振 点 的 方式 ,由 于 
过 绰 条 件 的 已 知 ,这 些 激 ( 测 ) 振 点 在 GFT 系统 中 实 
际 下 不 存在 的 。 
O 〇 CFT 系统 的 激 振 力 = Gf ,广义 力 为 
© Fp.=6(65,)=(6,.0,)"f,=0 7, 
式 钥 , 为 由 激 振 点 自由 度 表示 的 振 型 通过 插值 得 
GD 


AAA Z=0.0 
© 
CN 

虚拟 激 振 点 

> 

(a) 

CS 
a 

© 待 优化 激 振 点 


图 2 虚拟 激 / 测 振 点 技术 示意 图 


Fig.2 Virtual sensor-shaker technology view 


4 优化 算法 及 优化 流程 


选取 所 有 激 振 点 位 置 坐标 构成 优化 变量 x, 且 
X= fx zj, 对 于 梯形 机 辟 可 以 将 搜索 区 域 
影 为 正方 形 区 域 ( 见 图 2) ,构建 如 下 优化 问题 。 


定义 区 域 喘 射 时 ， 指 定 的 第 一 个 点 映射 
到 正方 形 区 域 的 (-1.1) 点 ， 剩 按 
顺 时 针 方 向 依次 映射 


义 力 和 原 系统 的 广义 力 相等 ,插值 振 型 $, 与 办 应 
该 尽量 接近 ,为 了 保证 插值 精度 , 则 需要 增加 激 振 点 
个 数 。 例 如 杆 面 模型 在 根部 固 支 ,因为 所 有 振 型 在 
根部 振 型 分 量 为 0, 如 果 在 根部 布置 激 振 点 则 会 减 
少 根部 区 域 的 插值 误差 ,大 大 提高 翼 面 振 型 插值 精 
度 。 同 时 ,由 于 在 根部 所 有 振 型 分 量 为 0( 即 固 支 没 
有 振动 ) ,因此 在 翼 面 根部 的 激 振 力 对 系统 的 广义 气 
动力 无 贡献 ,也 即 该 位 置 激 振 力 不 需要 物理 实现 。 
把 这 种 激 振 点 称 为 " 虚拟 激 振 点 ”, 即 这 些 激 振 点 只 
是 为 了 增加 插值 精度 ,而 实际 不 需要 施加 激 振 力 。 

对 于 己 面 根部 固 文 情况 ,同样 可 以 在 根部 增加 
“虚拟 测 振 点 ” ,这 些 测 振 点 在 实际 的 系统 中 不 必要 
存在 ,增加 的 目的 是 提高 振 型 的 插值 拟 合 精度 ,因此 
在 气动 力 计算 及 测 振 点 位 置 优化 中 是 有 必要 的 。 由 
于 玻 根 的 振动 为 0, 这 是 已 知 的 ,是 既 面 变形 的 有 效 
言 息 ,虚拟 测 振 点 可 以 利用 这 些 信 息 更 准确 地 完成 
恤 面 变形 插值 ,从 而 获得 更 加 准确 的 气动 力 。 示 意 
图 如 图 2 所 示 。 


F3 个 点 按 


_ a> 2 0 Coy 一 由 


min f,, = 了 (15) 


-1<x<1, C=1,2,.,n) 
采用 分 群 粒子 群 算法 进行 优化 搜索 。 相 比 标准 
型 粒子 群 算法 只 存在 一 个 固定 惯性 因子 的 搜索 群 


体 ,分 群 粒子 群 算法 将 粒子 群 分 成 知 干 个 粒子 子 群 ， 
每 个 子 群 对 应 不 同 的 惯性 因子 ,有 的 子 群 对 应 小 的 
惯性 因子 ,有 的 子 群 对 应 大 的 惯性 因子 ,这 样 就 对 各 
个 子 群 进行 了 分 工 , 惯 性 因子 小 的 子 群体 进行 局 部 
搜索 ,惯性 因子 大 的 子 群体 进行 全 局 搜索 。 这 样 做 
不 仅 保 证 了 整个 群体 全 局 寻 优 能 力 ,还 兼顾 了 群体 
的 局 部 搜寻 能 力 , 提 高 了 算法 的 寻 优 能 力 。 子 群 内 
的 粒子 在 设计 变量 空间 内 搜寻 运动 式 为 


k+l 
id 


7 = wo +ciri(P, — Xs) twry( Ps — Xx) 
(wg = Xi + Via 
式 中 :vi xi 分 别 为 第 迭代 步 第 i 个 粒子 优化 变量 
空间 内 第 d 维 速度 和 位 置 分 量 ;P; 为 第 i 个 粒子 个 
体 最 优 解 ;Py 所 有 粒子 的 全 局 最 优 解 ;1,7 为 0~1 
世间 的 随机 数 ;c 为 自我 信任 因子 ,2 为 种 群 信任 因 
了 一 般 取 ci =c =2;ow 为 惯性 因子 ,一 般 取 0.4 ~ 
0595 数 值 越 小 局 部 搜寻 能 力 越 强 | 。 

基于 广义 力 相等 的 GFT 系统 激 ( 测 ) 振 点 优化 

流程 如 下 (流程 图 见 图 3)。 


(16) 


原 芝 纹 剖 措 
分 相 


确定 激 “ 测 振 点 
日 标 振 弄 


沛 定 计 计 让 全 x 


图 3 激 / 测 振 点 位 置 优化 流程 
Fig.3 Optimization process of sensor-shaker location 
1) 首先 对 原 结构 进行 闸 振 分 析 , 获得 颤 振 分 文 
及 各 阶 模 态 参与 度 情况 。 
2) 确定 分 别 用 于 测 振 点 和 激 振 点 优化 的 目 
标 振 型 。 
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3) 确 定 适当 的 测 振 点 以 及 激 振 点 数量 ,并 给 定 
其 初始 位 置 。 

4) 分 别 以 测 振 点 位 置 和 激 振 点 位 置 作为 设计 变 
量 , 以 fiw = | $B -bb | 和 fi = 
6- 上 /9 作为 目标 函数 ,采用 优化 算法 
进行 优化 。 

5 ) 基于 激 振 点 和 测 振 点 的 优化 位 置 ,求解 GFT 
系统 的 控制 方程 ,对 GFT 系统 进行 颜 振 分 析 。 将 分 
析 结 果 与 原始 颤 振 系统 进行 比较 ,如果 精 度 满足 要 
求 则 结束 ;否则 ,返回 步 又 3) 。 


5 优化 算 例 与 试验 验证 


选取 一 个 梯形 平板 机 可 开 展 地 面 颤 振 试验 研 
究 , 对 其 激 ( 测 ) 振 点 位 置 进行 优化 ,机 辟 外 形 太 二 
如 图 4 所 示 ,材料 为 铝 合 金 。 


4 的 
[aa 


师 愉 :外 mm 


目 寓 图 


Bi HE 


到 4 梯形 平板 机 辟 外 形 
Fig.4 Trapezoid plate wing 
经 颤 振 分 析 ,结构 的 颤 振 机 理 主 要 为 前 两 阶 模 
态 耦合 ,因此 优化 过 程 中 选取 第 1、2 阶 振 型 作为 激 
振 点 和 测 振 点 的 目标 振 型 进行 优化 , 激 振 点 共 7 个 ， 
包含 3 个 虚拟 激 振 点 ,位 于 机 缀 根部 。 优 化 算法 采 
用 粒子 群 算法 ,优化 收敛 历程 如 图 5 所 示 。 可 以 看 
出 优化 目标 从 优化 前 的 0. 16 变 成 优化 后 的 0. 036， 
振 型 误差 降低 了 77.5% 。 优 化 前 后 激 ( 测 ) 振 点 位 
置 如 图 6 所 示 ,其 详细 位 置信 息 如 表 1 所 示 。 
对 优化 后 的 GFT 系统 进行 颜 振 分 析 ,得 到 的 广 
8 曲线 如 图 7 所 示 。 可 以 看 出 ,与 和 原 系统 对 比 , 优 
化 后 颤 振 分 支 和 颤 振 机 理 保持 不 变 。 优 化 前 后 GFT 
系统 颤 振 速度 和 频率 对 比 结果 如 表 2 所 示 。 优 化 后 


450 
颤 振 速 度 误差 减少 64% ,频率 误差 减少 70% 。 验 证 
了 优化 方法 的 有 效 性 。 
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6 优化 前 后 激 测 振 点 位 置 (4 个 点 情况 ) 


Fig.6 Comparison of the sensor-shaker location 


IV 


before and after optimization 
表 1 
Tab. 1 
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激 振 点 优化 前 的 位 置 坐标 
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tb) 优化 后 地 面额 振 系 统 广 g 曲 线 
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7 ”GFT 系统 优化 后 和 原始 闸 振 系统 V-g 图 
Fig.7 The V-g curve of the GFT-system 


and original flutter system 


激 / 测 振 点 位 置 优化 计算 结果 


The optimizing location for sensor-shaker 


激 振 点 优化 后 的 位 置 坐标 


振 点 
编号 % 村 1 bd * 各 
1 966.0 1 514.5 0.0 836.38 327. 63 0.0 
2 1 283.2 1 521.9 0.0 1 390.27 1 573. 64 0.0 
3 1 289.7 447.7 0.0 973.57 1 686.09 0.0 
4 530.9 445.8 0.0 660.31 1 143.55 0.0 
振 型 误差 0. 16 0.036 
优化 效果 振 型 误差 降低 了 [ (0. 16 -0.036)/0.16] x100% =77.5% 


表 2 优化 前 后 系统 颤 振 速度 和 频率 对 比 
Tab.2 Comparison of frequency and flutter 


speed before and after optimization 


GPT 系统 误差 减少 
参数 原 系统 “一 一 一 一 一 
优化 前 ”优化 后 /% 
颜 振 速 度 /(m .s-1) 118.56 130.63 114.24 64 
颜 振 频 率 /Hz 9.27 10.45 9.62 70 


基于 前 述 铝 合金 机 总 模型 ,利用 获得 的 优化 后 
激 / 测 振 点 位 置 , 搭建 了 地 面 闸 振 试验 系统 ,开展 地 
面 颤 振 试验 研 究 。 试 验 现场 如 图 8 所 示 。 实 时 气动 
力 解 算 模块 采用 频 域 气动 力 时 域 拟 合 方法 构建 时 域 
气动 力 模 型 ;采用 子 空间 辨识 方法 中 的 N4SID 算法 ， 
辨识 获取 系统 的 状态 空间 模型 进行 加 载 系 统 建 模 。 

试验 过 程 中 ,通过 不 断 调整 测试 来 流速 度 的 方 
式 进 行 试验 ,观察 结构 响应 收敛 和 发 散 情况 。 当 试 


验 速度 达到 127. 21 m/s 时 ,系统 响应 出 现 失 稳 临 界 
状态 ,如 图 9 所 示 ,对 应 振动 频率 为 8.45 Hz。 此 时 认 
为 系统 达到 闸 振 临界 状态 ,对 应 的 测试 速度 和 响应 的 
主 频 率 成 分 即 对 应 结构 的 颤 振 速 度 和 颜 振 频率 。 


| 


We 


@@ 激 光 测 振 仪 

图 8 GFT 试验 现场 
Fig.8 GFT test site 
1.5 
1.0 
& 0;5 


加 
CD 站 rs 人 
© 时 间 /s 时 间 /s 
CD (a) 临界 状态 (b) 发 散 状 态 
a 图 9 临界 失 稳 响 应 
OQ Fig.9 Critical instability response 


二 将 地 面 闫 振 结果 与 原始 闫 振 分 析 结果 对 比 , 斌 
验 的 颤 振 速 度 误差 为 7.3% , 赢 振 频率 误差 为 8. 8% ， 
EN 
地 葵 现 对 结构 的 颜 振 边界 的 测试 ,结果 的 误差 在 工 
程 接受 范围 ,证 明了 本 研究 提出 的 气动 力 模型 重 构 
方 力 的 正确 性 。 


6 结 论 


本 研究 平板 机 骂 GFT 系统 激 / 测 振 点 位 置 优化 计 
算 实例 验证 了 优化 方法 的 有 效 性 ,得 到 主要 结论 如 下 。 

1) 以 广义 力 相等 作为 优化 目标 考虑 了 原 系统 的 
颤 振 机 理 , 通 过 构建 加 权 系 数 ,能 够 保证 闸 振 分 支 和 
颤 振 机 理 不 发 生变 化 , 且 计 算 量 较 基于 闸 振 边界 等 
效 方法 小 。 

2) 虚拟 激 / 测 振 点 技术 能 够 提高 振 型 插值 精度 ， 
同时 有 效 减少 激 / 测 振 点 的 个 数 。 

3) 基 于 分 群 粒子 群 优化 算法 的 激 / 测 振 点 位 置 
优化 使 得 插值 点 位 置 分 布 更 合理 ,从 而 提高 GFT 系 
统 的 预测 精度 。 
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4) 基于 本 研究 方法 重 构 的 气动 力 构建 的 地 面 赢 
振 试 验 系统 ,开展 了 结构 颤 振 边界 测试 ,试验 结果 与 
分 析 结 果 吻 合 较 好 ,验证 了 方法 的 正确 性 。 
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